Structure of excited vortices with higher angular momentum in Bose-Einstein 

condensates 
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The structure of vortices in Bose-Einstein condensed atomic gases is studied taking into account 
many-body correlation effects. It is shown that for excited vortices the particle density in the 
vortex core increases as the angular momentum of the system increases. The core density can 
increase by several times with only a few percent change in the angular momentum. This result 
provides an explanation for the observations in which the measured angular momentum is higher 
than the estimation based on counting the number of vortices, and the visibility of the vortex 
cores is simultaneously reduced. The calculated density profiles for the excited vortices are in good 
agreement with experimental measurements. 
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When a superfluid is put into rotation, vorticity is split 
into discrete vortex lines rather than continuously dis- 
tributed as in the case of solid-body rotation Q, Q] . The 
dynamics, such as formation, reconnection and decay, 
of vortex lines is influenced by the microscopic struc- 
ture of the vortex core Q. Direct imaging of rotat- 
ing Bose-Einstein condensed atomic gases (BECs) has 
revealed, for the first time, the particle density profile 
of the vortex core in neutral superfiuids 0, 0, EIS Hi • 
The particle density is reduced in the core, and the den- 
sity dips in the measured images are used to identify 
the presence of vortex lines in the experiments. In the 
mean-field approximation, known as the Gross-Pitaevskii 
(GP) theory |3, |lfj , the predicted density at the center 
of the vortex core is exactly zero because the phase of 
the single-particle wave function at the vortex line is not 
well-defined. However, the observed density at the core 
center is finite and varies in a wide range with different 
experimental conditions. In many cases the core den- 
sities are quite large, and cannot simply be accounted 
for by the finite resolution of the imaging systems or 
by thermal excitations. Theoretically, finite densities 
at the core centers arise from quantum fluctuations of 
the locations of vortex lines due to many-body correla- 
tions HHIl [HQ- The fact that the core density is 
large and varies from measurements to measurements fur- 
ther suggests that these vortex lines may be in different 
vibrating modes ^ij rather than straight. The estimated 
density at the core center for a straight vortex line after 
taking into account quantum fluctuations is comparable 
to the density of particles depleted from the condensate 
in the ground state of a BEC m> 01 > which is not suf- 
ficient to explain the observed large core densities. In 
this Letter, the structure of excited vortices is studied 
using many-body wave functions that incorporate quan- 
tized motions of the vortex lines. It is proposed that the 
observed vortices with large core densities are in vibrat- 
ing modes corresponding to excited rotational states of 
BECs with an angular momentum higher than the sta- 
tionary GP states. This proposal is based on the results 



TABLE I: Ratios between the interatomic spacing, a, and the 
healing length, £, at the center of the traps in different ex- 
periments. In the Thomas-Fermi approximation, the ratio is 

given by a/£ = 2 (l5-K 5 Nm 3 alu)±U)\\/h 3 ) , where N is the 
total number of particles, m is the mass, a s is the s-wave scat- 
tering length, and u>± and are the trapping frequencies in 
the perpendicular and parallel directions to the rotation axis. 
The last column shows the ratio for liquid helium assuming 
that the vortex core radius is 1 A. 
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that the core density increases rapidly as the angular mo- 
mentum increases, and is supported by the experimental 
findings that the angular momentum per particle, aver- 
aged over an ensemble of one-vortex systems, is larger 
than h [lil Il7j . The existence of vortices with a higher 
angular momentum also provides an explanation for the 
observations in which the measured angular momentum 
is higher than the estimation based on counting the num- 
ber of vortices, and the visibility of the vortex cores is 
simultaneously reduced . The calculated density pro- 
files for the excited vortices are in good agreement with 
experimental measurements. 

The importance of many-body correlations in the vor- 
tex core can be seen in the following way. In the GP 
theory, the description of a vortex line is "classical" in 
the sense that the position and velocity of a vortex line 
can be simultaneously determined. A quantum descrip- 
tion would be required if the position uncertainty of a 
vortex line is comparable to the size of the vortex core. 
The position uncertainty of a vortex line is given by the 
interatomic spacing [bj, Il9| . because the total angular 
momentum of the system would be altered by one h if 
a vortex line fluctuates across an atom in the superfluid. 
The size of the vortex core is given by the healing length 
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which is determined through a balance between the ki- 
netic and potential energies for a density gradient. Al- 
though the healing length can be much larger than the 
interatomic spacing in the weakly interacting limit, they 
are often comparable with each other in realistic situ- 
ations (See Table __ for comparisons). As a result, the 
quantum fluctuations of vortex lines are generally impor- 
tant. We should note that it is possible to address the 
finite column densities within the framework of the GP 
theory by considering bent vortex lines Jljj. However, 
such considerations are limited to metastable states with 
broken rotational symmetry and with a lower angular 
momentum. 

For simplicity I consider a rotating BEC with only a 
single vortex point in two dimensions, but the formula- 
tion can be applied to a vortex line in a three-dimensional 
system. To obtain a quantum description for the vortex, 
many-body wave functions are constructed as linear com- 
binations of the GP wave functions parameterized by the 
location of the vortex [Til . The GP wave functions 
are used as basis states in analogous to the position space 
representation of a particle in quantum mechanics. For 
a A^-particle system, such a wave function is written as, 

N 

*(ri,---rjv) = / d 2 r F(r )n^Gp(ri;r ) , (1) 

i— 1 

where F(tq) is the weight function that represents the 
effective dynamics of the vortex, and V'Gp( r ; r o) is the 
normalized solution of the time-independent GP equa- 
tion with a vortex located at ro = (ro, Oq). When Eq. |QJ 
is generalized to describe a vortex line in three dimen- 
sions, the weight function becomes an effective "wave 
function" of a string rather than a particle. For the case 
of a straight vortex line with the zero-point motion, the 
wave function in Eq. JQ) is similar to the shadow wave 
functions used to study vortices in liquid 4 He |2l) . 

The weight functions for the energy eigenstates of the 
vortex can be solved by diagonalizing the many-body 
Hamiltonian within the sub Hilbert space spanned by 
the GP basis states. In the case that the system has 
a uniform ground state with a density p2D = c~ 2 , the 
weight functions are found to be 

F n ,i(r ) = Ar ntl x\ l \e- x ''*e ae ° 1 F 1 (a,b;a*), (2) 

where n > 0, n > I > — N, x — ^/nro/a, Af n j is the 
normalization constant determined by (^ n ,i\^n,i) = 1, 
a = -(2n - I - \l\)/2, b=\l\ + l and_i.Fi (a, 6; x 2 ) is the 
confluent hypergeometric function It is straightfor- 
ward to verify that the state \^ n ,i) is the eigenstate of the 
angular momentum with eigenvalue (iV + l)h, given that 
the phase of V'Gpforo) takes the form of 9 + f(9 — O ) 
in azimuthally symmetric systems. The weight func- 
tion .Fb,o( r o) is a Gaussian centered at the origin, and 
the state |^o,o) is the corresponding many-body state 



that includes the zero-point motion to the Hartree state 
of a centered vortex. The weight functions in Eq. (J2J 
have exactly the same form as the energy eigenfunctions 
of a charged particle in a constant magnetic field [23| . 
This is because the motion of a vortex is also driven 
by a velocity-dependent transverse force like the Lorentz 
force. For nonuniform systems such as BECs in harmonic 
traps, the exact solutions for the weight functions are not 
known. However, in the case that the radius of the sys- 
tem is much larger than the size of the vortex core and 
the vortex is near the trap center, the weight functions 
can be approximated by Eq. with /?2D equals to the 
local density of the ground state [14{ . Therefore, I will 
use the wave function in Eq. {_) with the weight func- 
tions given by Eq. (J_J to study the structure of vortex 
cores in trapped BECs. 

To use Eq. (JIJ, the first step is to set up the GP 
basis states. In order to solve the GP equation for a 
system with a given number of particles, it is conve- 
nient to choose the length unit to be the healing length, 
£ = ^/er 3 /87ra s , where a is the interatomic spacing of the 
ground state at the trap center, and a s is the s-wave scat- 
tering length characterizing the interaction strength be- 
tween particles. The ground state wave function tp ss (r) is 
first solved with the boundary condition ^ gs (Q) = 1. This 
particular choice of boundary condition and the length 
unit have the advantage that the normalization of the 
wave function is related to the total number of particles 
as 

N = _L J d 2 r|4 s (r)| 2 . (3) 

Then the GP state that has a centered vortex and the 
same number of particles is solved and written as 

^ GP (r;0) = g(r)j gs (r)e ie , (4) 

where g(r) is the amplitude ratio of the vortex state to 
the ground state. The normalized GP basis states with 
an off-centered vortex can now be approximated by 

V> G p(r;r ) « -^=-ff(|r - v \)^ gs (r)e^^ , (5) 
V No- 
where the phase satisfies the following set of linear dif- 
ferential equations, 

VxV<Xr;r ) = 2^(r - r ) , (6) 
V- [|V G p(r;r o )| 2 V0(r;ro)] =0. (7) 

The first equation specifies that the vortex is located at 
ro, and the second is the continuity equation which di- 
rectly follows from the GP equation. The general solu- 
tions of Eq. jnj can be written as (j)(r; r ) = $(r — r ) + 
<p(r;ro), where i?(r — ro) is the azimuthal angle of the 
vector r — ro, and 0(r;ro) is a single- valued function that 
remains to be determined by Eq. 0. 
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FIG. 1: (a) The integrated column density profile of a rotating 
BEC with a vortex at the center. The solid dots show the 
experimental data in Ref. [f|. The solid line shows the density 
profile of the excited state (n — 18) whose angular momentum 
per particle is about 1.05ft. I approximate the 3D column 
density using the 2D results as p co i(r) = p2o{r)%j) S s{r). The 
density profile of the GP state is shown as the dotted line for 
comparison. po(0) is the central density of the ground state, 
(b) The open circles show the scaling behavior of the central 
density with respect to the angular momentum. 



In the case that the vortex is close to the trap cen- 
ter, the correction term </>(r; ro) is only significant in the 
region, r > tq. In this region, the general solution of 
Eq. JHJl can be expanded as 



( r ; r o) 
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where the factor 1/j comes from the azimuthal angle 
#(r — ro). Since the density profile of the vortex state 
in this region is not far from the ground state density 
profile, we can replace |?/>Gp( r ; r o)| 2 in Eq. Q with the 
density profile of the ground state, po(r) — \ifj gs (r)\ 2 . By 
substituting Eq. (JSJ) into Eq. (JJJ , the differential equation 
for each cj>j(r) is obtained as follows, 



| P'o(r) d j 2 
r dr dr Pa{r) dr r 2 



ld_ d 



Mr) 



1 P'oir) 
r ]+1 po(r) 



, (9) 



where p' (r) = dpo(r)/dr. To solve Eq. @, two bound- 



ary conditions are needed. The boundary condition at a 
distance far away from the trap center can be obtained 
by noting that 
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(10) 



Thus, the asymptotic value of the first derivative of the 
solution is fixed by 
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(11) 



The other boundary condition can be chosen as <fij(0) = 
because the correction near the origin is small in a large 
system. With these two boundary conditions, Eq. @ 
is solved numerically to obtain <j)j(r), and thus the GP 
basis states. 

The second step is to find the overlap between two 
GP basis states because they are not orthogonal to each 
other. However, it is sufficient to calculate the overlap 
between two single-particle wave functions up to the sec- 
ond order in the separation between their vortex coordi- 
nates because the corresponding many-body wave func- 
tions are nearly orthogonal. The overlap between two 
single-particle wave functions is 
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(12) 



where cur and ai are numerical constants for a given 
system. The coefficient an generally scales logarithmi- 
cally with the system size, while the coefficient aj is not 
sensitive to the system size or local density variations, 
and is exactly given by 7r in an asymptotically uniform 
system |l4(. 

Finally, the density profile of the state \^f n i) is given 

by 



p n ,i(r) = N / d 2 r' / d 2 roK i (r[ 1 )F„j(ro)V'Gp(r;r )V'Gp(r;ro)exp 
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where the many-body overlap corresponding to the with the separation of the vortex coordinates in the large 
single-particle overlap in Eq. I|12|l falls off exponentially N limit. To compare this two-dimensional theory with 
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experimental data, we assume that the healing length £ 
and the interatomic spacing a at the trap center is the 
same as in the three-dimensional system. In the follow- 
ing, I will use the experimental conditions of Madison 
et al. 5] (The trapping frequencies are listed in Table 
P). For a BEC with N 3D = 1.4 x 10 5 particles, the heal- 
ing length is £ = 0.19 /im and the interatomic spacing is 
a = 0.17 /im. The corresponding number of particles in 
two dimensions given by Eq. J2J is about 350. The two 
coefficients in the overlap are clr — 5.84 and ai — 3.25. 

Although all states described by Eq. J2J) can exist, some 
of them are more stable and likely to be realized in ac- 
tual experiments than the others. A rotating BEC is typ- 
ically created by a stirring laser beam. The system is first 
driven by the laser stirrer for a period of time, and then 
relaxes to a stationary state. The laser stirrer breaks the 
rotational symmetry, and transfers both energy and an- 
gular momentum into the system. During the relaxation 
period, the energy can relax through collisional processes, 
but it is harder for the angular momentum to relax un- 
less the rotational symmetry is broken by other means. 
There is experimental evidence showing that the system 
can maintain its angular momentum for a long period 
of time ji^. This leads us to assume that the system 
will fall into the lowest energy state with the initially 
given angular momentum. As the angular momentum of 
the system increases, more and more vortices will form. 
However, it is energetically favorable to excite the exist- 
ing vortices if the change in angular momentum is small. 
The energy cost for introducing an additional vortex into 
the system is approximately tt ln(i?/£)?i 2 /m<7 2 , where R 
is the characteristic radius of the system, while the en- 
ergy cost for having an excited vortex is approximately 
proportional to nh 2 /ma 2 , where n is the angular mo- 
mentum increase and the proportionality constant is less 
than one 0, E| . Although adding another vortex would 
eventually become energetically favorable for larger val- 
ues of the angular momentum, we do not know what 
the critical value is because accurate comparison of the 
energies of different vortex configurations is complicated 
by the nonuniformity of the system and require further 
studies. 

In the case of a single excited vortex, the lowest en- 
ergy state for a given angular momentum is the \^ n ,n) 
state. Here the density profiles of the \^ n ,n) states are 
compared against experimental data. Figure ^a) shows 
that the data of Madison et al. 5] can be best fitted 
with the state that the angular momentum is increased 
by 5% (n = 18 compared to N = 350). The experimental 
data is shown by the solid dots. Only one centered vor- 
tex core is clearly visible, and the core density profile is 
significantly deviated from the density profile of the GP 
state (L z — Nh). Fig.^b) shows the trend that the core 
density increases with increasing angular momentum of 
the system (L z — (N + n)h). 

In summary, I have calculated the density profiles of 



different rotational states of a trapped Bose gas with one 
quantized vortex using many-body wave functions which 
are linear combinations of GP wave functions. The core 
density increases with increasing angular momentum of 
the system, which suggests that vortices with large core 
densities are in excited states with higher angular mo- 
mentum. The possibility of having excited vortices sug- 
gests that for given angular momentum, the lower energy 
rotational state can have number of vortices less than ex- 
pected in the mean-field theory. From the quantitative 
agreement on the density profile between the theoreti- 
cal prediction and the experimental measurements, it is 
suggested that the rotatingBose gas observed experi- 
mentally by Madison et al. [5j is in a state with angular 
momentum higher than one h per atom. 

The author thanks D. J. Thouless for valuable discus- 
sions, and K. W. Madison and J. Dalibard for providing 
data shown in Fig. ^ 
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